Nonlinear light localization around the core of a 'holey' fiber 
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We examine localized surface modes in the core of a photonic crystal fiber composed of a finite 
nonlinear (Kerr) hexagonal waveguide array with a central defect. Using a discrete approach, we 
find the fundamental surface mode and its stability window. We also examine an unstaggered, ring- 
shaped surface mode and find that it is always unstable, decaying to the single-site fundamental 
surface mode. A continuous model computation reveals that an initial vortex excitation (5* = 1) 
of small amplitude around the central hole can survive for a relatively long evolution distance. At 
high amplitudes, however, it decays to a ring configuration with no well-defined phase structure. 

PACS numbers: 42.81. Qb, 05.45.Yv, 42.65.Tg, 42.65.Wi 



INTRODUCTION 

Nonlinear propagation of fundamental Gaussian opti- 
cal beams has produced a rich variety of physical phe- 
nomena such as discrete and gap solitons in positive and 
negative periodic nonlinear media [IH3] . We can uncover 
an even wider range of novel nonlinear optical propaga- 
tion, by studying modes with different symmetries. One 
such mode is an optical vortex, which is an optical mode 
including a phase singularity at its centre. 

Optical vortices and their propagation have been stud- 
ied for their ability to trap and manipulate particles 0], 
and in the production of waveguides in atomic vapor [5] . 
The nonlinear propagation of vortex modes in the core 
of a photonic crystal fiber (PCF) have been studied the- 
oretically [tjHTU] , and there have been theoretical and ex- 
perimental study of vortex solitons in optically induced 
lattices [HHIl]. 

We take a new approach to the study of optical vor- 
tex propagation. Using an hexagonal array of nonlinear 
waveguides surrounding a solid core, we propagate an op- 
tical vortex in the waveguides adjacent to the core. We 
theoretically study the nonlinear propagation of vortex 
modes in this system, using both discrete and continu- 
ous models. 

Such structure is analogous to a liquid infiltrated pho- 
tonic crystal fiber (PCF) , which have been used to study 
nonlocal gap solitons [13] , the crossover from focusing to 
defocusing in a periodic array [15], as well as the possi- 
bility for selective infiltration for a range of interesting 
structures and applications [I6HT8] . 

By propagating a vortex mode in waveguides around 
the solid core of a PCF we can study vortex modes inter- 
acting with a surface, where the periodic structure meets 
a homogenous dielectric. Such states have been studied 
and observed in similar structures for single site excita- 
tion with Gaussian modes [THl HO] . 

This paper is organized as follows: Section II intro- 
duces the discrete model for an infiltrated PCF structure 



with an hexagonal geometry and a central defect making 
a solid core, in section III we introduce the continuous 
model of the same structure, focussing on the dynami- 
cal evolution of vortex excitations and finally, section IV 
concludes the paper. 

DISCRETE MODEL 

We consider a finite two-dimensional array of weakly- 
coupled nonlinear (Kerr) waveguides with hexagonal ge- 
ometry, with a missing waveguide at its center (Fig. 
[T] (a)). In the framework of coupled-modes theory, 
the electric field E(x, y, z) is presented as a superpo- 
sition of (single) transverse modes <fi(x, y) with ampli- 
tudes U(z) that vary slowly along the longitudinal di- 
rection: E(r,z) = J2Un(z)4>(r — ft), where f — (x,y) 
and ft — (n x ,n y ). These amplitudes obey the discrete 
nonlinear Schrodinger equation, 

i<N t +V H E4- l +7|E/n| 2 lfn=0 (1) 

where the sum is restricted to nearest-neighbors. The 
stationary solutions of Eq.Q have the form Uft(z) — 
Uffexp(i/3z), where Un obeys 

- /3 Uft + V + i\U H \ 2 U ri = 0. (2) 

We are interested in localized modes centered around 
the boundary of the solid core of the array. The simplest 
of these 'surface' modes is one centered on any of the 
six equivalent sites surrounding the missing guide. It is 
found by solving Eq. ^ using a direct extension of the 
Newton-Raphson method, starting from the decoupled 
(high- amplitude) limit, and performing a continuation 
process towards finite coupling values. For each mode 
found, we perform a standard linear stability analysis. 
Fig. flja) shows an example of a spatial profile for this 
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Figure 1. (Color online) Hexagonal waveguide array with cen- 
tral hole, (a) Example of a fundamental inner surface mode 
{P = 7.5, P = 8.33). Amount of shading denotes the distribu- 
tion of optical intensity, (b) Power vs. propagation constant 
curve for this kind of mode. Solid (dashed) curve denotes 
stable (unstable) portions. 



Figure 3. (Color online) Hexagonal waveguide array with 
central hole, (a) Example of an staggered 'ring' surface 
mode (/3 = 5.75, P = 47.95). Amount of shading denotes 
the distribution of the mode amplitude, going from white 
(amplitude=— 3) to dark gray (ampfitude=+3). (b) Power 
vs. propagation constant curve for this kind of mode. Solid 
(dashed) curve denotes stable (unstable) portions. 
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Figure 2. (Color online) Hexagonal waveguide array with cen- 
tral hole, (a) Example of an unstaggered 'ring' surface mode 
{P = 8.0, P = 37.4). Amount of shading denotes the distribu- 
tion of optical intensity, (b) Power vs. propagation constant 
curve for this kind of mode. This mode is always linearly 
unstable. 



uration, what are the decay channels for this mode? Will 
it transition to the low-power, single-site stable mode, or 
will it change into the staggered (stable) ring, or perhaps 
it will dissipate as radiation? To look for an answer, we 
follow the dynamical evolution of an initially completely 
localized ring mode configuration: Ua = Uq around the 
six sites surrounding the missing guide, Ua — other- 
wise. Long-time evolution of this mode over large prop- 
agation distance for a finite sample of N = 168 sites is 
shown in Fig. [4j Clearly, after a long transient behavior, 
where the diffracted beam bounces several times from 
the boundaries of the array, the beam becomes eventu- 
ally selftrapped in one of the six possible fundamental 
mode configurations. It is interesting to note that this 
selftrapping transition is quite abrupt, as evidenced in 
Fig. [5} 



kind of mode, along with its power content P = Wn\ 2 
vs propagation constant /? curve (Fig. [lib)). The curve 
obtained is typical of surface modes [21 and obeys the 
Vakhitov-Kolokolov stability criterion. In order to ap- 
proach something resembling a vortex-like mode, we con- 
sider next a higher-order mode, in the form of an unstag- 
gered 'ring' around the 'hole', with no phase difference 
(i.e., zero vorticity). An example of this high-power mode 
is shown in Fig. [2] along with its power vs. propagation 
constant curve. In this case, the mode is unstable for 
all values of its propagation constant. In fact, we find 
that most higher-order surface mode configurations are 
indeed unstable, with the exception of one: The stag- 
gered version of the ring mode (Fig. [3^ a)), where all am- 
plitudes around the hole are identical initially, but with 
a phase difference of n between nearest-neighbors around 
the ring. In this case, the mode is stable for initial am- 
plitudes (Fig. [3jb)) exceeding a given threshold. 

One interesting question at this points is: If we excite 
dynamically the unstaggered (i.e., unstable) ring config- 



CONTINUOUS MODEL 

Given the geometry of the array, it is conceivable that 
the addition of vorticity could stabilize this (unstable) 
ring mode. After all, we know that in two-dimensional 
square arrays, the addition of vorticity can stabilize some 
low-power modes that are otherwise, unstable [9]. 

In order to test this idea in conditions that are closer 
to an actual experiment, we simulate next the beam evo- 
lution in our structure by solving the continuous 2D non- 
linear Schrodinger equation for the slowly varying electric 
field envelope E: 

(3) 

Where = d\ + dy is the transverse Laplacian, 
D = X/Airn is the diffraction coefficient, 7 is the non- 
linear coefficient, A is the wavelength of light, n is the 
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Figure 4. (Color online) Evolution of unstable ring mode 
configuration over large propagation distance (a) 2 = 0, (b) 
z = 50, (c) z = 100, (d) 2 = 136, (e) z = 137, (f) z = 138(g) 
2 = 150, (h) 2 = 200. The initial amplitude in all six sites 
around the 'hole' is 2, with no phase differences. 



background refractive index, and G(x, y) is the refractive 
index profile defined numerically as a hexagonal lattice of 
circular holes with a diameter of d, pitch A which is the 
distance between the centre of two adjacent holes, and 
refractive index contrast An. While the general descrip- 
tion of the thermal nonlinearity is nonlocal [14] . here we 
use the simpler approximation of Kerr type nonlinearity. 

Such model is commonly used to investigate guidance 
properties in periodic arrays pQ. Using a hole diameter 
of d = 5/im and pitch A = 10/im, closely matching a 
commercial fibre (F-SM 15 by Newport) we are able to 
provide a theoretical basis for experimental observation. 

We propagate an input mode with profile E(z) = 
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Figure 5. (Color online) Evolution of unstable ring mode 
configuration over large propagation distance: Power content 
vs. longitudinal evolution distance, at the position of the 
center of the eventually self-trapped beam 




Figure 6. (Color online) (a) Linear propagation of an optical 
vortex mode with S — 1 at z — 2 cm (upper) and its phase 
(lower) in an hexagonal array of nonlinear waveguides. This 
linear mode diffracts as it propagates, (b) Nonlinear propa- 
gation of the same mode at z = 2 cm (upper) and its phase 
(lower) with 7 = 1 x 10~ 4 . The nonlinear mode propagates 
as a surface mode around a solid core defect. Pale dots on 
lower plots indicate waveguide location. 



rl s l exp(— ^2 + iS6), where E is the amplitude of the 
mode, S = 1 is the charge of the vortex, r, 9 and z are the 
cylindrical coordinates of the system, and w — y/2/\S\A 
is the width of the vortex mode. 

We find that even though the input mode is symmetric, 
linear diffraction causes some asymmetry in the output 
mode after z — 2 cm propagation (Fig. [6^a)) and index 
contrast An = 0.0032. The vortex phase is somewhat 
maintained at the output for linear propagation (one can 
pick a point in the cladding, and trace the phase in a 
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Figure 7. (Color online) (a) Linear propagation of an optical 
vortex mode with S = 3 at z — 2 cm (upper) and its phase 
(lower) in an hexagonal array of nonlinear waveguides. This 
linear mode diffracts as it propagates, (b) Nonlinear propa- 
gation of the same mode at z = 2 cm (upper) and its phase 
(lower) with 7 = 1 x 10~ 4 . The nonlinear mode propagates 
as a surface mode around a solid core defect. Pale dots on 
lower plots indicate waveguide location. 

circle around the core from to 2ir). This linear beams 
diffracts in the array as it propagates. Nonlinear output 
with 7 = 1 x 10~ 4 shows localization of the beam to the 
first ring of waveguides surrounding the solid core defect 
(Fig. |6|b)), and the loss of vorticity in the phase. Simi- 
lar to the discrete model, we see that the vortex mode is 
unstable when adjacent waveguides are not out of phase. 
These nonlinear modes are surface modes, largely con- 
fined to waveguides around the core for large propagation 
distance (z — 20 cm). 

To test the stability in the continuous model we next 
propagate a beam with 5 = 3, which satisfies the condi- 
tion of staggered phase between adjacent waveguides in 
the six waveguides surrounding the core. We see linear 
propagation is indeed more stable, and even more con- 
fined in the form of a ring mode, even though the vortex 
phase is lost (Fig. |7^a)). The linear beam diffracts as it 
propagates in a similar fashion to S = 1 modes. In the 
nonlinear regime we see the mode begins to break up as 
the staggered phase is lost (Fig. |7[b)). Again these non- 
linear modes are surface modes, although this time more 
confined to waveguides around the core over large propa- 
gation distance (z — 20 cm), due to the initial staggered 
phase profile. 

The linear and nonlinear modes produced with this 
5 = 3 input mode are more strongly confined to the 
waveguides adjacent to the core defect, when compared 
to modes produced with an input with 5=1. While 
the 7T staggered phase between all six waveguides around 
the core is lost in the nonlinear propagation, the sites 
which have the highest intensity maintain this staggered 



phase. The initial vorticity of the 5 = 3 beam seems to 
stabilize the linear output into a vortex which survives 
for long propagation distances. In the nonlinear regime 
the vorticity is lost, but the ring mode structure is main- 
tained with a somewhat staggered phase between some 
waveguides. 



CONCLUSIONS 

In summary, we have examined the localized surface 
modes around the core defect of a PCF surrounded by a 
hexagonal array of nonlinear waveguides. We find that 
the stable modes in both the discrete and continuous 
models have a staggered phase profile for the six waveg- 
uides surrounding the core. Ring shaped surface modes 
are studied in the discrete model and shown to always 
decay to a single site fundamental surface mode. The 
continuous model shows a similar decay of the surface 
modes and loss of vorticity in the phase at high nonlin- 
earity. 

It is suggested that this work could be performed in 
an experimental setting using a liquid infiltrated PCF. 
One needs to be careful in choosing the parameters of 
the fibre, and the nonlocal character of the nonlinearity 
must be taken into account |14j . 
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